Symmetric and antisymmetric exchange anisotropies in quasi-one-dimensional 

CuSe 2 5 as revealed by ESR 
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We present an electron spin resonance (ESR) study of single-crystalline spin chain-system 
CuSe205 in the frequency range between 9 GHz and 430 GHz. In a wide temperature range above 
the Neel temperature TV = 17 K we observe strong and anisotropic frequency dependence of a 
resonance linewidth. Although sizeable interchain interaction Jic ~ 0.1J (J is the intrachain inter- 
action) is present in this system, the ESR results agree well with the Oshikawa-Afneck theory for 
one-dimensional 5=1/2 Heisenberg antiferromagnet. This theory is used to extract the anisotropics 
present in CuSe20s. We find that the symmetric anisotropic exchange J c = (0.04 ± 0.01) J and 
the antisymmetric Dzyaloshinskii-Moriya (DM) interaction D — (0.05 ± 0.01) J are very similar in 
size in this system. Staggered-field susceptibility induced by the presence of the DM interaction is 
witnessed in the macroscopic susceptibility anisotropy. 

PACS numbers: 75.10.Pq, 75.30.Et, 75.30.Gw, 76.30.Fc 



I. INTRODUCTION 

Magnetism of quasi-one-dimensional (ID) S = 1/2 sys- 
tems is often well described by an isotropic Heisenberg 
Hamiltonian. Low dimensionality enhances quantum 
fluctuations, which in turn suppress long-range order- 
ing. Small interchain interactions present in real systems 
usually stabilize long-range magnetic ordering at low, 
but finite, temperatures. The ground state of such sys- 
tems is, however, very sensitive to the presence of small 
anisotropy of exchange interactions between spins, frus- 
tration and/or defects. Depending on a local symmetry, 
both symmetric and antisymmetric anisotropic exchange, 
i.e. the Dzyaloshinskii-Moriya (DM) interaction, 1 can be 
present. In ID systems where staggered g tensor and/or 
staggered DM interaction is present, an applied mag- 
netic field induces a staggered field, which opens a gap in 
the excitation spectrum. 2,3 Determination of the leading 
anisotropic terms of the spin Hamiltonian thus represents 
an important milestone in understanding these materials. 

Magnetic resonance techniques are a very powerful tool 
for addressing the above points and the type of a ground 
state in such systems. 4,5 One of the most appropriate and 
sensitive methods for determining the presence of small 
anisotropies is electron spin resonance (ESR), 6,7 since 
anisotropic exchange interactions broaden the otherwise 
exchange narrowed ESR line. 8 Kubo-Tomita (KT) the- 
ory is a well-established method of the linewidth analysis, 
however, it is unfortunately limited only to high tempera- 
tures T ^> J. 9 Analysis of the ESR linewidth has evolved 
considerably in the last 15 years, especially regarding 
the S = 1/2 ID Heisenberg antiferromagnet (HAF). 10 " 12 
In case of a staggered DM interaction, its contribution 
to the ESR linewidth is of the same order of magni- 



tude as the symmetric-anisotropy contribution at high 
temperatures, 12 despite the fact that latter is expected 
to be smaller, being a higher order perturbation correc- 
tion to exchange coupling. 1 The perturbation theory cal- 
culations of the ESR line in S = 1/2 ID HAF have been 
extended to the entire temperature range only for the 
case of symmetric anisotropic exchange. 13 Relatively re- 
cently Oshikawa and Affleck (OA) employed field-theory 
methods to derive a general low-temperature ESR re- 
sponse for half-integer spin ID HAF. 10,11 Their predic- 
tions for the ESR linewidths were experimentally verified 
in a system with the dominant DM interaction. 14 Sepa- 
rately, low temperature theories 10,11,13 were recently suc- 
cessfully applied to systems with symmetric anisotropic 
exchange. 15 However, quantitative analysis of the low- 
temperature ESR linewidths for the realistic systems 
with both anisotropies present has so far been limited. 

Among the most studied ID systems arc copper oxides 
in which magnetism originates from Cu 2+ spins 5* = 1/2 
and superexchange is usually mediated through Cu-O-Cu 
pathways. Copper(II) diselenium(IV) pentoxide is a new 
ID copper oxide. It crystallizes in the monoclinic space 
group C*2/c. 16,17 The structure consist of S = 1/2 Cu 2 + 
chains running along the crystallographic c axis (Fig. 1). 
Each Se 4+ ion carries one lone pair of electrons which 
plays a role of "chemical scissors". 17 The most relevant 
exchange interactions should be those shown in figures 
1(a) and 1(b). The dominant intrachain interaction J 
is mediated through the double Cu-O-Se-O-Cu super- 
superexchange (SSE) paths. From a phenomcnological 
(structural) point of view, one may expect significantly 
smaller interchain (IC) coupling Jjc than J. The ar- 
rangement of the CUO4 plaquettes and the Cu-O-Se-O- 
Cu bridges is such that it prevents sizeable next-nearest 



2 




y A 



although, the disagreement below T max remains notice- 
able, especially for the chain direction. 18 Despite be- 
ing predominantly ID, the system orders magnetically 
at Tjv = 17 K 18 due to sizable interchain interac- 
tions, as evidenced by Jic ~ TV- The susceptibility 
anisotropy below Tjv is consistent with an antifcrromag- 
netic type of spin arrangement, with spins oriented per- 
pendicular to spin chains. 18 Raman scattering measure- 
ments showed that spin-spin correlations emerge below 
ps 110 K, coinciding with T max , and that the system 
is dominated by enhanced classical spin dynamics as a 
consequence of a rather strong interchain interaction. 20 
So far CuSe205 was treated as isotropic Heisenberg spin 
system. 18,20 In this work we concentrate on neglected 
magnetic anisotropics and obtain the anisotropic spin 
Hamiltonian of the system by simultaneously modelling 
the angular, the temperature and the frequency depen- 
dence of the ESR linewidth with the OA theory. With 
the obtained anisotropy terms we significantly improve 
the agreement between measured and modelled suscepti- 
bility anisotropy below T = J. The results presented 
here show that both the symmetric and the antisym- 
metric anisotropic exchange need to be taken into ac- 
count in case of CuSe2 05. In ID cases when both types 
of anisotropies are present and comparable, temperature 
and frequency dependent ESR measurements are invalu- 
able for their quantitative assessment. 

II. SPIN HAMILTONIAN OF CuSe 2 5 



FIG. 1. (color online). Linear-chain crystal structure in 
CuSeaOs. (a) Orange plaquettes represent Cu04 rectangles 
and grey spheres represent Se atoms, (b) Intrachain super- 
superexchange path with exchange coupling J (solid lines) 
and dominant interchain path with exchange coupling Jic 
(dotted lines), (c) Local coordinates of the staggered g ten- 
sor, the staggered DM vector D = {D a *, 0, D c ) and the axis 
of the symmetric anisotropic exchange J c . Dotted line repre- 
sents the chain. 



neighbour coupling along the chain. There are two short 
IC exchange paths, but one, with the modest size of 
Jjq ps 0.1 J, is expected to be significantly stronger than 
the other [dotted lines in Fig. 1(b)]. 18 

The ID magnetic character of the system is reflected 
in bulk magnetic susceptibility displaying a broad max- 
imum at T max = 100 K. 18,19 Above T max it can be well 
modelled already in the simplest ID Heisenberg limit, 
yielding the dominant intrachain exchange J = 157 K. 18 
This value is in-line with the Curie- Weiss temperature 
®cw — 165 K and was also successfully theoretically 
accounted for by density-functional-theory (DFT) cal- 
culations, predicting J = 165 K and Jic = 20 K. 18 
Taking into account the interchain couplings from the 
DFT calculations, quantum Monte Carlo (QMC) cal- 
culations could further improve the agreement between 
experimental susceptibility and theoretical predictions, 



The spin Hamiltonian, which describes the spin-spin 
interactions in a quasi- ID spin system in the applied mag- 
netic field H, is 

H = Ui S o + Uae + H D M + Uz, (1) 

where 

Hiso = J ^ S; • Sj+i + Jic Si ■ Sj, (la) 

i <i,j> 

Hae = Jn ^ ] S i S i+1 , (lb) 

i 

K M/ =^D l .(S ! xS l+1 ), (lc) 

i 

Hz = -Hb S * ' & ' H ( ld ) 

i 

Hi S o is the isotropic Heisenberg interaction with intra- 
chain exchange J and interchain exchange Jic- The first 
sum in H.i 80 runs over spins along the chain and the 
second over spins along the chain S, and over nearest- 
neighbouring spins Sj from two neighbouring chains. %ae 
is the symmetric anisotropic exchange with symmetry 
axis n and relative magnitude S = J n /J, Hdai is the 
antisymmetric DM anisotropic exchange term 1 with the 
site-dependent DM vector [see Fig. 1(c)]. Hz is the 
Zeeman term where g; is the g tensor for site i. 
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The orientation of the anisotropic-exchangc symmetry 
axis n may be nontrivial to determine when two or more 
incquivalent sites are present. 21 In CuSe20s the situation 
is further complicated by the fact that inequivalent sites 
A and B [Fig. 1(c)] are bridged by the complicated Cu- 
O-Se-O-Cu SSE path. Nevertheless, below we show that 
the anisotropic exchange symmetry axis coincides with 
the chain direction, i.e. n = c in Eq. (lb) [sec Fig. 1(c)]. 

The general form of the DM vector D = (D a * , 0, D c ) 
is imposed by a two-fold rotational axis along b passing 
through the middle of each intrachain Cu-Cu bond. 1 In 
addition, the symmetry of CuSe2C>5 is such that the DM 
interaction is staggered, i.e. ±D [Fig. 1(c)]. 

The crystal symmetry also dictates the g tensor to be 
staggered for A and B sites [Fig. 1(c)]. For these two 
sites we thus split gi in Eq. (Id) into uniform, g„, and 
staggered, g s , component and from now on use 



TABLE I. Results of linear fit AH = A + B T of high tem- 
perature linewidth in 400 — 550 K range. 



Ss ' — Su ± gs ■ 



III. EXPERIMENTAL 



(2) 



Single crystalline CuSe20s samples were synthe- 
sized by standard chemical vapour transport method, 
as described previously, and characterized by X-ray 
diffraction. 17 

The ESR experiments were performed at X-band 
(9.4 GHz) and at high-frequencies (HF) between 50 GHz 
and 430 GHz on single-crystalline samples. Temperature 
dependence between 4 K and 550 K was measured with 
temperature stability better than ±0.05 K. X-band mea- 
surements were performed on a home-made spectrometer 
equipped with the Varian TEM104 dual cavity and the 
Oxford Instruments ESR900 cryostat. High-frequency 
ESR was performed using custom-made transmission 
type spectrometers at NHMFL facility at Tallahassee, 
Florida. 22 In all cases the ESR spectra were fitted to a 
single Lorentzian line. 

Magnetic susceptibility anisotropy was determined 
from torque magnctomctry measurements performed on 
a home-made torque apparatus. The resolution of the 
magnetometer is better than 10~ 4 dyn cm. Measure- 
ments were performed in magnetic field of 8 kOe in the 
temperature range 2 - 330 K. 



IV. RESULTS 
A. X-band ESR 

The room-temperature (RT) angular dependence of 
the X-band ESR spectra (inset of Fig. 3) is shown in 
Fig. 2. It reveals g-factor values for the three crystallo- 
graphically relevant directions, g a * = 2.064, gi, = 2.140 
and g c = 2.226. Taking into account the Cu site sym- 
metry we determine the principal eigenvalues of the g- 



Direction 


a* 


b 


c 


A (Oe) 


220 ± 20 


223 ± 7 


313 ±15 


B (Oe/K) 


0.70 ±0.05 


0.61 ±0.01 


1.10 ±0.06 



tensor g x = 2.064, g y = 2.089 and g z = 2.277, which are 
considerably higher than those previously used in theo- 
retical calculations. 18 The principal axes of the g ten- 
sor with respect to the CUO4 plaquette are shown in 
Fig. 1(c). As expected from the crystal structure, lo- 
cal crystal-field symmetry at the copper site is close to 
being uniaxial with the local anisotropy axis pointing in 
the direction perpendicular to the CUO4 plaquette, i.e., 
being tilted by a = 32° from the c axis around the a* 
axis. We take into account that the total measured g 
tensor is in the strongly exchanged narrowing limit given 
by g = (g A ± g B )/2. The uniform g u and the staggered 
component g s defined in Eq. (2) thus have the following 
form in the a* b c frame 

'2.064 
2.140 
2.226y 

'0 
0.084 
^0 0.084 

The measured g factors are temperature independent in 
the paramagnetic state and only slightly increase in the 
vicinity of Tjy ■ 

At RT, the ESR linewidth anisotropy is pronounced in 
the a*c plane and marginal in the plane perpendicular 
to the chains (a* 6 plane). At lower and higher temper- 
atures, however, the anisotropy of AH in the a* b plane 
slightly increases, as can be seen in Fig. 3. The tem- 
perature dependence of the ESR linewidth along crys- 
tallographic directions a*, b and c was measured from 
T N = 17 K to 550 K (Fig. 3). Above w 200 K the 
linewidth increases linearly with temperature for all ori- 
entations. From the value of J it is not expected that 
spin-spin correlations would persist up to 550 K, so we 
attribute this linear dependence to the phonon assisted 
spin-lattice broadening. 23 We fit our high-T data to the 
phenomenological expression AH(T) = A ± B T where 
AH ph = B T is the phonon- induced line broadening. 
Parameter A is the temperature independent exchange- 
narrowed linewidth, as predicted by the KT theory in the 
T > J regime. 9 The parameters obtained from the fits 
are summarized in Table I. The values of parameters A 
and B slightly depend on the temperature range of the 
fit, which has been accounted for in the parameter errors 
listed in Table I. 

Subtracting the phonon contribution AH ph from the 
raw data AH results in the corrected linewidth AH° = 
AH — AH ph , which is determined by the spin-spin in- 
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FIG. 2. (color online), (upper panel) Angular dependence of 
the g-factor at RT in X-band. (lower panel) Angular depen- 
dence of the linewidth at RT i n X-band. Solid lines are fits 
to g = JgJ, sin 2 9 + g\ h cos 2 6. 



teractions only. At high temperatures AH C is given by 
the temperature independent coefficients A (Table I) and 
starts to gradually decrease with decreasing temperature 
below RT. At w 100 K there is a crossover to even 
steeper decrease with decreasing T, which is a conse- 
quence of the evolution of spin-spin correlations below 
Tmax- 18 At T « 25 K the linewidth has a minimum and 
then starts to sharply increase with decreasing tempera- 
ture. Below T/v ESR spectrum disappears at X-band. In 
principle, the increase of the linewidth in the vicinity of 
the Neel point is expected due to critical slowing down 
of spin fluctuations. 24 In this case the theory predicts 
that the linewidth becomes sensitive to the resonance 
frequency, 25 which was experimentally seen as a decrease 
of the linewidth with increasing frequency. 26 This is in 
contradiction with measurements in CuSe2 05, as can be 
seen in Fig. 4. Here the observed linewidth increases with 
frequency even at temperatures well above T/v, thus sug- 
gesting some other origin of field dependent linewidth. 
OA showed that in ID S = 1/2 HAF the linewidth in- 
creases with decreasing temperature for T <C J, if stag- 
gered fields are present. 10,11 Since in CuSe20s both the 
staggered g tensor and the staggered DM interaction are 
potential sources of staggered fields, we now turn to high- 
frequency ESR results. Nevertheless we note that large 
increase in X-band linewidths below « 22 K is reflecting 
the critical fluctuations in the vicinity of T/v . 



B. High frequency ESR 

Temperature dependences of the ESR linewidth mea- 
sured at 240 GHz for crystal directions a*, b and c are 
shown and compared to the X-band data in insets of 
Fig. 4. For comparison we also show linewidth measured 
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FIG. 3. (color online). Temperature dependence of X-band 
ESR linewidth measured along three crystallographic direc- 
tions. Solid lines represent linear fits to high temperature 
data (see text). Inset: RT spectra for a* and c direction. 

at 112 GHz for direction b. In contrast to the X-band 
data, at higher frequencies the ESR signal is observable 
even below T/v = 17 K, showing a clear anomaly in the 
linewidth at the transition temperature. 

It is immediately clear that the magnetic field strongly 
affects the ESR linewidths along all crystallographic di- 
rections. By far the largest effect is seen for the mag- 
netic field along the c direction. To quantitatively anal- 
yse the temperature and frequency (field) dependence of 
the linewidth we first subtracted the high temperature 
phonon contribution AH ph determined at X-band fre- 
quencies. In doing so, we exploit the fact that this con- 
tribution is field independent. The resulting linewidths 
AiT, which will be discussed below, are shown in the 
main panels of Fig. 4. 

The temperature region where the OA theory is strictly 
applicable is T/v « T C J, although experiments have 
been successfully analysed also for T < J. 14 In CuSe20s 
J « 160 K and the system orders antiferromagnetically 
at T/v = 17 K, so we limit our analysis to the tempera- 
ture range 22 K < T < 100 K. Increasing the lower limit 
by a few Kelvins does not influence the extracted param- 
eters. For the Hamiltonian given by Eq. (1) the theory 
for the ESR linewidth in case of symmetric anisotropic 
exchange predicts 11,13 

A^) = ^4r, (4) 

where e = 2 applies when magnetic field is along the 
anisotropy n-axis and e = 1 otherwise. Low-temperature 
logarithmic correction can be neglected in the investi- 
gated temperature range. 11 Eq. (4) predicts linear in- 
crease of the ESR linewidth with T, which is exactly 
what we observe for AH C for all three directions and 
T > 25 K (Fig. 4). Thus we conclude that the symmet- 
ric anisotropic exchange is present in CuSe2 05. 

Clear deviations from the linear temperature depen- 
dence of AH C in the low-temperature region imply that 
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FIG. 4. (color online). Temperature dependence of the cor- 
rected ESR linewidth AH C measured at 9.4 GHz, 112 GHz 
and 240 GHz for magnetic field orientations along all three 
crystallographic directions. Solid lines are fits to Eq. (6). In- 
sets: temperature dependence of the raw-data linewidth mea- 
sured at different frequencies. 



the staggered field is also present. Its contribution to the 
linewidth is given by 1 ' n 



AH sf {H, T) = 0.69 g /jb 



ks J 



'In 



(5) 



where the staggered field h s = c s H is proportional to 
the applied field H and the anisotropic staggered field 
coefficient c s . 

Since there are no cross terms between the symmetric 
and the antisymmetric anisotropic exchange, the total 
ESR linewidth for ID system described by Hamiltonian 
(1) is given by 

AH (H, T) = AH + AH ae (T) + AH sf (H, T). (6) 

AHq is the temperature independent linewidth in the 
high-T limit, which - according to Refs. [10, 11] - con- 
tains both the field-independent and the field-dependent 
contributions. 

Simultaneous fitting of both the temperature (Fig. 4) 
and the frequency dependence (Fig. 5) of AH C to the 
model given by Eq. (6) allows us to obtain parame- 
ters AHq, S and staggered field coefficients c s .i for all 
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FIG. 5. (color online). Frequency dependence of the cor- 
rected ESR linewidth AH C at T = 25 K. Solid lines are fits 
to Eq. (Al). 



three magnetic field orientations. The details of the 
more involved analysis of the frequency dependence are 
given in Appendix A. Both linewidth dependencies can 
be fitted (Figs. 4 and 5) with the single set of pa- 
rameters: c s>a , = 0.009 ± 0.004, c Sjfc = 0.021 ± 0.004, 
c s . c = 0.062 ± 0.006, S = 0.04 ± 0.01. Simulations also 
yield that the c axis is the symmetric-anisotropy axis, 
namely n = c in the Hamiltonian (lb). The worse agree- 
ment is found for X-band data below rj 25 K, which we 
attribute to the vicinity of the phase transition overshad- 
owing the effect of the induced staggered field. To affirm 
the correctness of the above fitting parameters we now 
turn to the low-temperature angular dependence of the 
linewidth at high frequencies. 



V. STAGGERED FIELD ANALYSIS AND THE 
DM VECTOR 

In CuSe2 05 the crystal structure implies that both 
the staggered g tensor and the staggered DM vector are 
present, giving the staggered field 



1 
2J 



D x g, 



H. 



(7) 



The staggered field contribution to the ESR linewidth is 
determined by the magnitude of h s , which is proportional 
to the staggered field coefficients c s . For the cases when 
magnetic field is aligned along the three crystallographic 
axes, in CuSe2C>5 the staggered field coefficients have the 
following form 



■*s,a* 



\Dc\g a * 
2g b 



( D c g b \ 2 { D a ,g b + 2g bc ^ 



/ V 2 5a* J 2gc 

Zgbc - D a ,g c 



(8) 



256 



6 





1 1 1 1 1 1 1 1 


1 1 1 

Q 


i i i | i i | i 

□ data 






IraP T 


\ D * ' 




/ »V 








it y 




;.\ --" D - 




/ w 








/T / 












\?\ 

\ "A 








\ \ 












w 




\ & 







. i , 


i . 



-90 -60 -30 30 60 90 

0(deg) 

FIG. 6. (color online). Angular dependence of ESR linewidth 
measured at 240 GHz and 25 K in the o*c plane. Solid and 
dotted lines are calculated linewidths using Eq. (6) and D + = 
(-0.044, 0, +0.0255) J and D~ = (-0.044, 0, -0.0255) J, 
respectively. In both cases 5 = 0.04. See text for details. 



The general expression for c s for arbitrary direction of 
magnetic field can be found in Appendix B. 

Using the experimentally determined g tensor [Eq. (3)] 
and values of c s ^ 1 we can now obtain also the DM vector 
from Eqns. (8). Solving c s a * and c s j, for D gives four 
possible solutions D± = (-0.114 ± 0.015, 0, ±0.024 ± 
0.014) J and D± = (-0.043±0.015, 0, ±0.024±0.014) J. 
Alternatively, solving c StC and c Si b gives two real values: 
Df = (-0.045 ±0.011, 0, ±0.027 ±0.014) J. Since DM 
vectors and Dg nearly coincide, we take this solution 
as the correct DM vector. 

Finally, we note that the average of T>f and D^, i.e the 
DM vector D = (-0.044 ± 0.010, 0, ±0.0255 ± 0.010) J, 
accounts also for the angular dependence of the linewidth 
measured at 240 GHz and 25 K (Fig. 6). Unfortunately, 
the present data does not allow us to determine also the 
sign of the c-component of D as calculated linewidth for 
both cases describe data equally well. 



VI. DISCUSSION AND CONCLUSIONS 

The main experimental finding of this work is the de- 
termination of anisotropies in the ID HAF CuSe205. 
Both anisotropies, the symmetric anisotropic exchange 
and the antisymmetric DM interaction are sizeable and 
comparable. The former anisotropy is expected to open 
a gap in the excitation spectrum already in zero mag- 
netic field. This gap is expected to decrease with increas- 
ing field. 27 However, in contrast at high magnetic fields 
the staggered field contribution will dominate and define 
the energy gap. In this limit, the excitation gap is ex- 
pected to scale with magnetic field, similarly as in other 
prototypical ID HAF systems, like Cu-benzoate 2,3 ' 28 or 
Cu-pyrimidine dinitrate. 29,30 The crossover between the 
two regimes has to be addressed in the future from both 
theoretical as well as experimental point of view, since it 
may lead to unconventional static and dynamic magnetic 



properties. 

Presence of both anisotropies should explain the unre- 
solved issue of poor fitting of magnetic susceptibility be- 
low T max . ls Staggered magnetic anisotropy is responsible 
for the additional anisotropy in the magnetic susceptibil- 
ity for T < J. For this reason we measured the tempera- 
ture dependence of the magnetic susceptibility anisotropy 
&Xba* = Xb ~ Xa* and Ax C a* = Xc - Xa* by means of 
torque magnetometry (Fig. 7). The measured magnetic 
susceptibility anisotropies clearly cannot be explained 
within simple ID HAF Bonner-Fisher susceptibility 31 ' 32 
(Fig. 7). In calculations we used J = 161 K, while in- 
terchain interactions with Jjc = 0.1J and number of 
nearest neighbouring chains z = 2 have been taken into 
account in the mean-field approximation. 32 The quality 
of the fit significantly improves when staggered suscepti- 
bility Xs,i(T) is taken into account. Xs,i(T) is expressed 
as 2,3 

where i represents the direction of the applied mag- 
netic field. Using above given staggered susceptibility 
we expand a quantitative agreement with susceptibility 
anisotropy data down to 50 K with somewhat increased 
staggered field parameters, c s , a * = 0.024, c s ,b = 0.047 
and c SjC = 0.13, the discrepancy between the staggered 
field parameters obtained from ESR and from susceptibil- 
ity being similar to what was observed for Cu benzoate 3 . 
Thus determined magnetic anisotropies even manage to 
describe the plateau in Axca* below T max as well as 
mimic the sudden upturn in A%6a* below 50 K. 

In conclusion, we studied quasi-one-dimensional 
CuSe205 with temperature, frequency and angular de- 
pendent ESR. ESR linewidth analysis within the OA 
theory for S = 1/2 ID HAF allowed us to obtain 
the values of symmetric anisotropic exchange interaction 
J c = (0.04 ± 0.01) J and antisymmetric DM interaction 
|D| = (0.05±0.01) J. CuSe205 thus appears to be an ex- 
tremely interesting system where both anisotropies are of 
similar strength, which has some profound effects on the 
ground state and magnetic susceptibility anisotropies. 
Present results thus challenge detailed investigations of 
the low-energy excitation spectrum and the staggered 
spin susceptibilities by inelastic neutron diffraction or lo- 
cal probe nuclear magnetic resonance techniques. Finally 
we stress, that this work also demonstrates how to sys- 
tematically approach ESR data when both anisotropies 
are present in the ID system. 

Appendix A: Frequency dependence of the linewidth 

at T = 25 K 

Frequency dependence of the linewidth measured at 
25 K for all three crystallographic directions is shown in 
Fig. 5. To fit the data to the Eq. (6), we need to re- 
solve Eq. (6) into frequency-dependent and frequency- 
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FIG. 7. (color online). Susceptibility anisotropy measured 
in the a* b plane (upper panel) and a*c plane (lower panel). 
Lines represent comparison between susceptibility anisotropy 
obtained from ID HAF with interchain interactions (dotted 
lines) and when staggered susceptibility is included (solid 
lines). See text for details. The error bar is shown on the 
right of the data. Vertical line represents the ordering tem- 
perature Tn- 



independent contributions. The frequency dependent 
contributions are expected to depend on v 2 , so we fit 
our data for given direction i to: 

AHi(T, v) = Oi(T) + ft(T) v 2 , i = a\ b, c, (Al) 



where 



a;(T) = AH ,i(v -- 
ft(T) v 2 = AH ,i(v) 



0) + AH ae>i (T), (Ala) 
AIL,, ,\u. I ;. (Alb) 



The temperature independent part AHq in Eq. (6) is also 
split in two parts, the frequency independent AHoj(i> = 
0) and the frequency dependent AHq^(u). Our intention 
is to obtain the parameters 8 and c s . However, both pa- 
rameters a and ft obtained from the fit consist of two 
unknown parts. The difficulty arises from the fact that 
both AHq^{u = 0) and AHq^(u) need to be determined 
first. The origin of these contributions was discussed by 
OA. 11 From their analysis AHo t i(u = 0) is expected to 
be negligibly small because it is proportional to D 4 /J 3 . 
If we neglect it we obtain 8 = 0.03 which is in good agree- 
ment with 8 = 0.04 ± 0.01 obtained from fits of tempera- 
ture dependence of the linewidth. To determine the fre- 
quency dependent high-T term AHq^(v) we assume that 
the 9 GHz data is the zero-frequency linewidth limit and 
subtract it from the 240 GHz linewidth to obtain the 
frequency-dependent part, AH$(v = 240 GHz). From 
Eq. (Alb) we see that by subtracting this from ft (25 K) ■ 



TABLE II. Results of the fit of the frequency dependence of 
the linewidth at T = 25 K shown in Fig. 5 to Eq. (Al). 





a 





AH (240 GHz) 


Cs 




(Oc) 


(10" 3 Oc/GHz 2 ) 


(Oe) 




a* 


60 ± 20 


0.7 ±0.2 


20 ± 5 


0.010 ±0.002 


b 


120 ± 18 


2.9 ±0.2 


20 ± 5 


0.026 ± 0.001 


c 


62 ±61 


17.3 ±0.8 


90 ±5 


0.068 ± 0.002 



(240 GHz) 2 we obtain AH sf (v = 240 GHz, 25 K). Once 
AH sf (v = 240 GHz, 25 K) is obtained the staggered field 
coefficients can be calculated by using Eq. (5). The re- 
sults are summarized in Table II and correspond, within 
errors, to the ones obtained from the fits of the temper- 
ature dependence of the linewidth. 



Appendix B: Angular dependence of the linewidth 

at T = 25 K and v = 240 GHz 

Assuming the most general staggered DM vector al- 
lowed by symmetry, D = (D a * ,0, D c ) J, the following 
expression is obtained for the staggered field coefficient 
c s from Eq. (7) 



D c gt, sin 9 sin tp 
2^ 

/ (2g bc - Da* g c ) cos 9 + D c ga' cos ip sin t 
\ 2g b 

( (Da* gb + 2ffb c ) sin 9 sin tp 
{ 2g~ c 



(Bl) 



g tensor components can be directly read from Eq. (3), 9 
is the polar and tp the azimuthal angle (we take x = a*, 
y = b and z = c). Eq. (Bl) gives the coefficient c s for 
some general direction of the applied field. Experimen- 
tally determined staggered field coefficients (Sec. IV B) 
are equated with (Bl) in the following way: c S;a * = 
c s (9 = n/2,tp = 0), c s . b = c s (9 = n/2,p = tt/2) and 
Cs.c = Cs{Q = 0), which gives result (8) quoted in the 
main text. 

Measured angular dependence of the linewidth shown 
in Fig. 6 can be compared to Eq. (6) without any free 
parameters using expression (Bl) for the staggered field 
parameter by setting tp = 0. We assume the following 
angular dependence of parameter e in (4): 6 = 1 + cos 2 9 
which gives e = 1 for a* direction and e = 2 for c, in ac- 
cord with the previously obtained results. Temperature 
independent correction AHq in Eq. (6) is also angular 
dependent. We simulate this dependence in the follow- 
ing way: AH (6) = AH^ a , + (Aff ,c - AJTo.oO co ^ e ^ 
where AH , a * = 40 Oe and AH 0iC = 450 Oe are the 
temperature independent corrections obtained from fits 
of Eq. (6) to the temperature dependence of linewidth 
along the a* and the c axis, respectively. We also 
take 8 = 0.04 obtained from the same fits, and D ± = 
(-0.044, 0, ±0.0255) J. 
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